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ABSTRACT 

We search for signs of rotation in the subsystems of the Milky Way and M31 that are de- 
fined by their satellite galaxies, their globular cluster populations, and their BHB stars. A 
set of simple distribution functions are introduced to describe anisotropic and rotating stellar 
populations embedded in dark haloes of approximate Navarro-Frenk- White form. The BHB 
stars in the Milky Way halo exhibit a dichotomy between a prograde rotating, comparatively 
metal-rich component ([Fe/H] > —2) and a retrograde rotating, comparatively metal-poor 
([Fe/H] < — 2) component. The prograde metal -rich population may be associated with the 
accretion of a massive satellite (~ 10 9 M Q ). The metal-poor population may characterise the 
primordial stellar halo and the net retrograde rotation could then reflect an underestimate in 
our adopted local standard of rest circular velocity Go. If 0o is ~ 240 kms -1 then the metal- 
poor component has no rotation and there is a net prograde rotation signal of w 45 kms" 1 
in the metal-rich component. There is reasonable evidence that the Milky Way globular clus- 
ter and satellite galaxy systems are rotating with (v^,) ~ 50 kms -1 and {v$) « 40 kms" 1 
respectively. Furthermore, a stronger signal is found for the satellite galaxies when the an- 
gular momentum vector of the satellites is inclined with respect to the normal of the disc. 
The dwarf spheroidal satellites of M3 1 exhibit prograde rotation relative to the M3 1 disc with 
(v<f,) w 40 kms -1 . We postulate that this group of dwarf spheroidals may share a common 
origin. We also find strong evidence for systemic rotation in the globular clusters of M31 
particularly for the most metal-rich. 

Key words: galaxies: general - galaxies: haloes - galaxies: kinematics and dynamics - galax- 
ies: individual: M31 - dark matter 



1 INTRODUCTION 

Our own Milky Way Galaxy and its nearest neighbour Andromeda 
(M31) provide testing grounds for theories of galaxy formation and 
evolution. Studies of local stellar populations explore the complex 
processes which built up the galaxies. The orbits of tracer popula- 
tions bear the imprint of the infall and subsequent accretion history 
of the dark matter halo. In a monolithic collapse model (see e.g., 
lEggen et alJI 1962b . the halo and the disc of the galaxy are drawn 
from the same population. Any rotational signature in the outer 
halo is thus aligned with the angular momentum of the disc. How- 
ever, in a hierarchical picture, substructure arriving late at large 
radii may have little connection to the disc. Its evolution is fur- 
ther complicated by processes such as tidal disruption, disk shock- 
ing and dynamical friction. For example, prograde rotating sys- 
tems are more susceptible to dynamical friction than retrograde 
dOuinn & GoodmarJll986h . iNorris & Rvanl Jl989h suggested that 
the subsequent disintegration of these accreted fragments by tidal 
forces may create a net retrograde asymmetry in the remnant stellar 
halo population. The presence, or absence, of a rotating signal can 
provide key insights into the formation history of the halo. 
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There have been a number of previous kinematic studies 
of halo populations in the Milky Way Galaxy. For example, 
lFrenk&Whitehl980l) found that the globular cluster system has 
a velocity distributio n that is isot ropic about a net prograde rotation 

extended this work by restricting 
attention to halo globular clusters with [Fe/H] < — 1 an d found 
a sligh tly weaker net prograde motion (v rot ~ 50 kms^^. lNorrisI 
q 19860 used a sample of metal-poor objects ([Fe/H] < —1.2) and 
found evidence for weak prograde rotation, although he found no 
kinematic difference between the clusters and halo stars, suggesting 
these are part of the same population. 

The isotropic velocity distribution of the globular clus- 
ters contrasts with the highly anisotropic dis t ributions usually 
found for Population II star s dWoollevI 1 19781 : lHartwickl 1 19831: 
iRatnatunga & Freemanl Il985l) . Cosmological simulations gener- 
ally predict increasi ng radial anisotron y with Galactocentric radii 
for halo stars (e.g. lAbadi et"al] 120061) . as a consequence of the 
accretion of infalling structure. Observati onally, the situation i s 
less clear-cut with radial l y bias ed (e.g. Ichiba & Beersl |2000|) . 
isotropic (e.g. ISirko et al.l | 2004bl) and even tangentially biased 
dSommer-Larsen et al .11 199m velocity distributions claimed for the 
outer stellar halo. 

More recently, a number of authors have advocated a dual 
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halo structure, which suggests that the stellar halo was formed by 
at least two distinct phases of accr etion events. The evidence in- 
clude s studies of the spati al profiles dHartwickll987l : lKinman et"aD 
1 1994 | Miceli et"al] 120081) . as well as kinematic studies indicat- 
ing net retrograde rotation in the outer parts ( :Majewski| 1 1 9921 ; 
ICarnev et al.ll 19961 IWilhelm et al Jfl996l: iKinman et al.ll2007l) . For 



example, Carollo et alj 1 2007 ) claimed a dual halo structure con- 
sisting of a weakly prograde rotating inner halo and a retrograde 
rotating outer halo of lower metallici ty ([Fe/H] < —2). Ho wever, 
this idea is inconsistent with w ork by Chiba & Beerj d2000l) on lo- 
cal neighbourhood stars and bv lSirkoetal.N2004bh on distant BHB 
stars who both find no evidence for rotation in the outer halo. 

This somewhat confused picture is partly caused by the diffi- 
culties of the task. Full kinematic analyses of tracer populations are 
hampered by small sample sizes and lack of full information on the 
phase space coordinates. Often, work has been confined to com- 
paratively local samples, which may not be unbiased tracers of the 
overall population. This, though, is beginning to change with more 
rece nt studies on BHB stars reaching ou t to distances of ~ 80 kpc 
(e.g. lSirko et alj|2004bl : lxue et alj|2008r) . albeit with only the line- 
of-sight velocity available. In fact, several surveys have the promise 
to provide much larger samples of halo tracers with well-defined 
kinematical parameters. This inc ludes ongoing p rojects like the 
Sloan Digital Sky Survey (SPSS, lYork et al .1120001) , and the Sloan 
Extension for Galactic Understanding and Exploration ( SEGUE, 
iNewberg & Sloan Digital Sky Survey Collaboratior|2003l) . as well 

il\2~ 



; future projects like the Gaia satellite dTuron et al.t2005n. and the 
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Large Synoptic Survey Telescope (LSST, 

There has also been a lively debate in recent years as to 
the possibility that the satellite galaxies of the Milky Way lie 
in a rotation a lly-su pported disc. The idea may be traced back to 
iLvnden-Belll d 19831) . who first suggested the Magellanic Clouds, 
Ursa Minor and Draco may have been torn from a single pro- 
genitor, a gigantic gas-ric h proto-Mage l lanic C loud. Working with 
more complete datasets, iKroupa et alj d2005l) claimed that satel- 
lite galaxies occupy a highly inclined disc and that this is at odds 
with the predictions of cosmological simulations. The available 
proper motion measurements have been used t o constrain the an- 
gular momenta o rientations of the satellites (e.g. lPalma et al.l2002l ; 
iMetz et al1l2008l) . indicating evidence for some sort of coherent 
motion, possibly rotation. Even if true, it remains disputatious as to 
whether this can be reconciled within the cold dark matter frame- 
work of structure formation dMetz et all 1200 8l ; iLibeskind et alj 
l2009h . 

Rotational properties of halo populations can therefore pro- 
vide clues as to their origin and evolution, and may allow us to iden- 
tify associations sharing a common formation history. With this in 
mind, we develop some simple and flexible distribution functions 
that can model rotating stellar populations embedded in dark ha- 
los. Cosmological arguments su ggest that dark halos have a univer- 
sal Navarro-Frenk- White form dNavarro et alj|l99q) . whilst many 
stellar populations in the outer halo are well-approximated by a 
power-law. Section 2 provides a flexible set of simple distribution 
functions that are powers of energy E and angular momentum L 
and which can be adapted to include rotation. We apply the mod- 
els to the satellite galaxy, globular cluster and BHB populations of 
the Milky Way (Section 3) and M31 (Section 5), with the implica- 
tions of our results for the Milky Way given particular attention in 
Section 4. Finally, Section 6 sums up. 



2 DISTRIBUTION FUNCTIONS 

2.1 The Even Part of the Distribution Function 

A collisionless system can be described by a phase-space distribu- 
tion function (henceforth DF), F. The probability that a star occu- 
pies the phase space volume d 3 a: d 3 v is given by F(x, v)d 3 x d 3 v. 
DFs are a valuable tool for studying steady-state systems as they re- 
place the impracticality of following individual orbits with a phase- 
space probability density function. The part of the DF that is even 
in v<p fixes the density, whilst the part that is odd in fixes the 
rotational properties. 

To construct the DFs, we assume a steady-state spherical po- 
tential and a density profile for our halo populations. Of course, 
these are tracer populations moving in an external gravitational 
field generated by the dark matter halo rather than the self- 
consistent density generated through Poisson's equation. For sim- 
plicity, we use simple power-law profiles for the density and poten- 
tial, namely p oc r~ a and $ oc r 1 , where a and 7 are constants. 

The density profiles of dark matter halos in cosmological 
simulations resembl e the Navarro-Frenk- White or NFW profile 
dNavarro etai]|l996l) . Exterior to r ~ 20 kpc, we can approxi- 
mate the NFW potential by a power-law profile, as shown in Fig.Q] 
Here, we have used one of the NFW models applicable to the Milky 
Way from lKTypin et all d2002h . and have normalised th e power-law 
profile to reproduce the local escape speed found by ISmith etHll 
d2007l) . We see that a power-law with index 7 = 0.5 is a good 
approximation at large radii (r > 20 kpc). At small radii, the ap- 
proximation breaks down, but in this regime, the bulge and disk 
components become important. In any case, our main application 
is to tracer populations that reside far out in the halo. The dark halo 
mass within a given radius is 



M(< r) = 



1/2 



1/2 



2G 



(1) 



where <3>o is the normalisation at the solar radius Rq . For a virial ra- 
dius of 250 kpc, the halo mass enclosed is ~ 8 x lO n M0, in good 
agreement with estimates f or the Milky Way Galaxy and M31 (e.g . 
Evans & Wilkinsonl |2000| : IXue et alj I2OO8I : iBatl 
Watkins et alfeOld) . 



attaglia et al] |2005t 



Previous authors have adopted power-law den sity profiles fo r 
halo popula tions with 2 < a < 5. Fo r exam ple, lHarrisI dl976h . 
IZinnl d 19851) and iDiorgovski & Mevlanl d 19941) find a power-law 
with a = 3.5 is a good fit to the Milky Way globular cluster pop- 
ulations. Classical studies have fo und a = 3.5 fo r stellar popula- 
tions in the Milk y Way halo (e.g. lFreemanlll987l) . More recently, 
iBell etal.ld2008l) find values of a in the range 2 — 4, although they 
caution against the use of a single power-law due to the abundant 
substructure in the stellar halo. The density profile of the system 
of Milky Way satellite galaxies is poorly kno wn due to the small 
sampl e size and incompleteness at large radii. lEvans & Wilkinsonl 
< 200lh find a = 3.4 is a good approximation for r > 20 kpc whilst 
I Watkins et al. adopt a power law of a = 2.6. Similar pro- 

files h ave been used for the M31 halo populations. ICrampton et all 
dl985h argue that the radial profile of the M31 globular cluster 
system is similar to that of the Milky Way and so it is rea son- 
able to adopt t he sam e power- law. lEvans & WilkinsorJ d2000l) and 
I Watkins et al. use a values of 3.5 and 2.1 respectively, for 

the M31 satellites. As a convenient summary of all this work, we 
will use a density power-law index of a = 3.5, but we discuss the 
effects of varying a in Section[3] 

Armed with these simple forms for the potential and density, 
we can give the velocity distribution in terms of the binding energy 
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Figure 1. V esc as a function of radius. The power law approximation devi- 
ates from a NFW profile at r < 2 kpc. The NFW pa rameters applicable 
to the Milky Way were taken from lKlypin et alj j2002h . The normalisation 
of the power law is chosen to match the local escape velocity derived in 
ISmithetalJi2007l) . 



Figure 2. The mean streaming velocity (iu ) as a function of Galactocen- 
tric distance for 7 = 0.5 and a = 3.5. Solid and dotted lines correspond 
to maximum prograde and retrograde rotation respectively. The profile flat- 
tens at large radii. Blue, red and green lines represent isotropic, radial and 
tangential velocity distributions respectively. 



E = <&(r) — h(vf + v\ + vf os ) and the total angular momentum 
L = h\ + Ly + E\ as 

F cvcn (E,L)<xL~ 2l3 f(E) 

where 



(2) 



g(-y-2) . a 3 



(3) 



Here, /3 is the Binney anisotropy parameter dBinnev & Tremainel 
1 1987b . namely 



= 1 



(v 2 o) + (vl) 



(4) 



which is constant for the DF s of the form of eqn {2j. These DFs are 
discussed in greater detail in lEvans et al.l J 1997h . 



2.2 The Odd Part of the Distribution Function 

So far, our DFs describe non-rotating populations embedded in 
dark halos. We now devise an odd part to the distribution function 
which generates a one-parameter family of rotating models 



F odd = (1 - r/)tanh(L*/A)F e , 



(5) 



where r\ is a constant. The case r\ = describes maximum pro- 
grade rotation, whilst 77 = 2 describes maximum retrograde rota- 
tion. Here, A is a 'smoothing' parameter to ease numerical calcu- 
lations and soften the Heaviside function. 

To illustrate the effects of the parameter 77, we compute the 
mean streaming velocity (v<p), which is analytic for 7 = 0.5 



_ 2?r(f-fl)r(2a- 4/3 + 1) 
W 7rr(l - P)T{2a - 4/3 + §) 



(6) 



Table U| lists (v^,) values at r — 50, 100 kpc for three different 
anisotropy parameters. Fig.[2]shows the radial profile of the stream- 
ing motion for maximum prograde (solid line) and maximum ret- 
rograde (dotted line) distributions. 



V 


P 


W, 50 k P c)kms 1 


("0,100 k P c)kms 1 


0.0 


(-0.5, 0, 0.5) 


(82, 72, 53) 


(69,61,45) 


0.5 


(-0.5, 0, 0.5) 


(41,36, 27) 


(34, 30, 22) 


1.0 


(-0.5, 0, 0.5) 


(0, 0, 0) 


(0, 0, 0) 


1.5 


(-0.5, 0, 0.5) 


(-41,-36, -27) 


(-34, -30, -22) 


2.0 


(-0.5, 0, 0.5) 


(-82, -72, -53) 


(-69, -61, -45) 



Table 1. Typical values evaluated at r = 50 kpc and r = 100 kpc 
as a function of the rotating parameter r] and the anisotropy parameter /3, 
[7 = 0.5, a = 3.5] 



2.3 Properties 

It is rare that we possess full six-dimensional phase space informa- 
tion for any star. However, we can still use a DF by marginalising 
over the unknown components. For example, consider a case where 
the spatial position and line of sight velocity are well known, but 
the proper motions are uncertain. In this case, we marginalise over 
vi and Vb to obtain the line of sight velocity distribution (LOSVD): 

F(l, b, d, vios) — J J F(l,b,d,vi,v b ,v los )dvidv b . (7) 

Once 7 is fixed as 0.5, there remain the free parameters, /3, a and 
V- 

Some properties of the distribution functions are shown in 
Fig.[3] The top left panel shows F (I, b, d, V\ os ) as a function of v\ os 
at three different distances from the Sun (for the case 01 — 3.5). 
An arbitrary line of sight toward the Leo constellation is chosen for 
illustration, (I, b) = (230°, 50°). Closer to the Galactic Centre the 
escape velocity is larger, so the velocity distributions are broader. 
Tangential orbits (/3 < 0) have narrower radial velocity and broader 
tangential velocity distributions. The reverse is true for radial orbits 
03 > 0). 

For distant objects in the Milky Way, v\ os ~ v r and so the 
line of sight velocity distribution reflects the radial distribution. At 
small radii, the trend is not so simple and depends on the particular 
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Figure 3. Top left: The dependence of the distribution function on ?Ji os at three different distances in the non-rotating case. A random line of sight is chosen 
(toward the Leo constellation) and the density power law is set to a = 3.5. Blue, red and green lines represent isotropic, radial and tangential velocity 
distributions respectively. Top right: The distribution function for three different density profile power laws in the non-rotating case. The distance is fixed at 
d = 80 kpc. Middle left: The line of sight velocity distribution when rotation is included. A prograde rotating case is shown (r) = 0.1) for two different 
longitudes. At large distances little difference is seen from the non-rotating case. Middle right: As the adjacent panel, but now the distance is fixed to 
d = 30 kpc. More pronounced rotation is seen at smaller distances. Bottom left: The line of sight distribution function for the prograde case at three different 
distances. The same line of sight as the top panels is chosen. Very small deviations from a non-rotating distribution are seen. The tangential and radial velocity 
distributions are shown by green and red lines respectively. Bottom right: The longitudinal velocity, Vi, distribution with the same parameters as the adjacent 
panel. Rotation is much more apparent in this case and the distribution is almost entirely skewed to negative u; . 
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line of sight in question. The top right panel of Fig. [3] illustrates 
the dependence of F(l, b, d, v\ os ) on the power law index of the 
density profile, a. In a given gravitational potential, more extended 
tracer populations (smaller a) exhibit larger velocity dispersions 
than more centrally-concentrated ones. 

The middle panels of Fig. [3] show the changes to the line 
of sight velocity distributions caused by the introduction of rota- 
tion. At large distances (d ~ 80 kpc) the rotating distribution is 
hardly distinguishable from the non-rotating case (as v\ os fa v r ). 
As expected, radially biased velocity distributions show the small- 
est deviation from the non-rotating case, but even a tangentially bi- 
ased system may only show small differences. At smaller distances 
(d ~ 30 kpc) the profile depends on the particular line of sight as 
illustrated in the middle right panel. 

The bottom panels show the line of sight velocity distribution 
and the longitudinal velocity distribution for prograde rotating sys- 
tems (rj — 0.1). We see that with just the line of sight velocity, the 
rotation of a system may be difficult to constrain. In particular, if 
the line of sight component of velocity is largely due to the radial 
component, then we can only poorly constrain rotation. However, 
applications where the line of sight velocity has significant contri- 
butions from the tangential velocity components, as for example in 
the case of M31, can give a more robust prediction for the rota- 
tion. The bottom right panel shows the distribution as a function 
of vi. In this case, the distributions show much more pronounced 
differences from the non-rotating case. For this particular line of 
sight and with rj = 0.1, the distribution is almost entirely skewed 
to negative vi . This has important implications for the relevance of 
accurate proper motion measurements. The components of motion 
perpendicular to the line of sight can give conclusive evidence for 
rotation, whilst the deductions from the line of sight velocities rely 
heavily on the particular line of sight. 



3 APPLICATIONS: THE MILKY WAY GALAXY 

In this section, we apply our distribution functions to halo popu- 
lations in the Milky Way Galaxy. Observed heliocentric velocities 
are converted to Galactocentric ones by assuming a circular speed 
of 220 kms -1 at the position of the sun (Ro = 8.5 kpc) with a so- 
lar peculiar motion_£C7 J J/jy}=(ll.l, 12.24, 7.25) kins" 1 (recently 
updated bv lSchonrich et alj|20ich . Here, U is directed toward the 
Galactic centre, V is positive in the direction of Galactic rotation 
and W is positive towards the North Galactic Pole. 

We set the power law indices for the potential and density to 
7 = 0.5 and a = 3.5 respectively. We allow the anisotropy f) 
and rotation rj parameters to be estimated by maximum likelihood 
analysis of the data. For example, when we have full distance in- 
formation and the line of sight velocity, we construct the likelihood 
function from the LOSVD: 

N 

L(r ] ,j3)=Y t fo&F(li,b i ,di,v loai ,r,,l3), (8) 
i=i 

where N is the number of objects in the population. Eqn gives 
the two dimensional likelihood as a function of /3 and rj. As we are 
mainly interested in constraining rotation, we sometimes fix /3 and 
find the likelihood as a function of rj. 

3.1 Milky Way BHB stars 

Blue Horizontal Branch (BHB) stars are excellent tracers of halo 
dynamics as they are luminous and have an almost constant abso- 



lute magnitude (within a certain colour range). We construct a sam- 
ple of BHB stars from the Sloan Digital Sky Survey (SDSS) DR7 
release. To select BHBs, we impose constraints on the (extinction 
corrected) u — g and g — r colours, the surface gravity (g 3 ) and the 
effective temperature (T c g), namely: 

0.8 < u - g < 1.4 
-0.4 < q - r < 

(9) 

2 < log( 5s ) < 4 
7250 < T cff /K < 9700 

These cuts are similar to those used by previous authors (e.g. 
lYannv et al.ll200(]|.ISirko et alj|2004ij) and are implemented to min- 
imise contamination by blue stragglers and main sequence stars. 
We only consider stars wit h \z\ > 4 kpc, as we wish to exclude 
disc stars. IXue et al.l J2008 1 ) compiled a sample of 2558 BHB stars 
from SDSS DR6 using Balmer line profiles to select suitable candi- 
dates. Our larger sample (as we use the DR7 release) of 5525 stars 
includes « 88% of the Xue sample. 

Heliocentric distances are evaluated photometrically by as- 
suming BHBs have an absolute magnitude of M g = 0.7 in the 
g band Q. The absolute magnitudes of BH B stars are slightly af- 
fected by temperature and metallicity (e.g. ISirko et alj[2"004ah but 
we expect that this variation makes a negligible difference to our 
distance e s timate s. We compare our distance estimates to those of 
IXue et al. I fe008h . who take into account such variations,and find 
agreement to within 10%. 

We remove contamination of stars belonging to the Sagittarius 
dwarf galaxy by masking out tracers lying in the region of the lead- 
ing and trailing arms. Our final sample contains 3549 BHB stars 
(with r > 10 kpc) and so provides a statistically representative 
population to study the outer stellar halo. The heliocentric velocity 
errors are « 5 kms" 1 on average, at larger radii the error increases 
due to the decreasing brightness of the source. The maximum error 
in the velocity is ~ 20 kms -1 . 

W e use the metallicities derived from the IWilhelm et al.l 
d 19991) analysis (the 'WBG method) as this method was adopted 
specifically for BHB stars. The method uses the Call K lines as 
a metal abundance indicator and is applicable in the range — 3 < 
[Fe/H] < 0. Only a small fraction (N ~ 300) of stars do not 
have assigned metallicities or are outside the applicable metallicity 
range. 

Fig. [4] shows that the overall population (purple contour) 
has slight net retrograde rotation (77 = 1.1) and has an al- 
most isotropic distribution. By splitting the sample into metal- 
poor ([Fe/H] < —2, blue contours) and 'metal-rich' ([Fe/H] > 
—2, red contours) subsamples, we find retrograde and pro- 
grade signals for the two subsets respectively. The mean stream- 
ing motions for the metal-poor and metal-rich samples are 
= -(35 ± 10)(r/10kpc)- 1/4 kms" 1 and (i> ) = (21 ± 
ll)(r/10 kpc) -1 / 4 kms -1 respectively. We remove the radial de- 
pendence by averaging over the volume interval (10 < r/kpc < 
50) to find (ity) = -(25 ±7) kms -1 and («^) = (15 ±8) kms -1 
for the metal-poor and m etal-rich populations respectively. These 
results can be compared to lCarollo et al .1 d2007h who find a rotating 
retrograde outer halo and a weakly prograde rotating inner halo. 
Our analysis only applies to the outer halo, but instead of a sin- 
gle population with net retrograde rotation, we find evidence for 

We check that changing our assumed absolute magnitude from M g = 
0.7 to M g = 0.55 makes a negligible difference to our dynamical analysis. 
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Figure 4. lcr confidence regions for Milky Way halo BHB stars in the 77, /3 plane from a maximum likelihood analysis. Arrows indicate prograde, retrograde, 
tangential and radial distributions respectively. The overall population (N = 3549 with r > 10 kpc) is shown by the purple filled contour. Metal poor 
([Fe/H] < —2, N = 1135) and metal 'rich' ([Fe/H] > —2, N = 2125) subsamples are shown by blue and red contours respectively. The paler shades 
show the radial bin 25 < r/kpc < 50 whilst the darker shades show the radial bin 10 < r/kpc < 25. 



two separate popu lations split by meta llicity (at [Fe/H] ~ —2). 
In agreement with ICarollo e t al. (2007|) (who focus on halo stars 
within the solar circle) we find evidence for a net retrograde rotat- 
ing metal-poor halo component, but we find that the outer halo also 
contains a more metal-rich slightly prograde-rotating component. 

We can make a rough estimate of the luminosity associated 
with these metal-rich and metal-poor remnants assuming 40% of 
the metal-poor sample (77 ~ 1.4, iV P oor = 1135) and 20% of 
the metal-rich sample (77 ~ 0.8, iVrich = 2125) are rotating. 
The fraction of rotating stars can be scaled to the overall popula- 
tion by calculating the total number of BHBs in a spherical vol- 
ume. This requires a normalisation factor for our density profile 
(p oc r~ ). We adopt the normal isation applicable to RR Lyrae 
stars found bv lWatkins et al .1 d2009l) and assume the ratio of BHBs 
to RR Lyrae is approxima tely 2:1. F inally , using photometry for 
globular clusters given in lAn et al. I d2008h . we find the relation 
Nbhb/L ~ 10 _3 Lg 1 (where Nbhb is found by applying a cut 
in colour-colour space - see eqn|5J. We find for both the metal-poor 
and metal-rich rotating populations a luminosity L ~ 10 8 L@. For 
a mass to light ratio of M/L ~ 1 (appropriate for luminous satel- 
lites - see Fig. 9 in lMateolll998l ) this gives a total mass estimate 
of M ~ 1O 9 M . This rough calculation suggests that these ro- 
tating signals could be remnants of two separate accretion events 
caused by massive satellites with different metallicities and orbital 
orientations. 

The different shaded contours in Fig. [4] show different radial 
bins. The paler shades are for the range 25 < r/kpc < 50 whilst 
the darker shades are for the range 10 < r/kpc < 25. Splitting the 
sample into radial bins shows the rotation signal is present at all 
radii (albeit with less confidence at greater distances due to smaller 



: 10 < r/kpc < 25 

2 ~ 25 < r/kpc < 50 



— 




-2 7 i - 

2.0 2.5 3.0 3.5 4.0 4.5 5.0 
a 

Figure 5. The variation in the anisotropy parameter ft with the density pro- 
file index a for Milky Way halo BHB stars. The solid line and dashed 
lines are for the radial bins, 10 < r/kpc < 25 (N = 2312) and 
25 < r/kpc < 50 (N = 1127) respectively. More steeply declining 
density profiles (larger a) give more radially anisotropic /3 values. Arrows 
represent tangential and radial distributions respectively. 
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Figure 6. Line of sight velocity dispersion as a function of Galactocentric 
radii for Milky Way halo BHB stars. The metal poor ([Fe/H] < -2) and 
metal rich ([Fc/H] > —2) components are given by the blue and red lines 
respectively. An obvious cold feature is apparent in the metal rich popu- 
lation at r ~ 25 kpc. The dotted line is the average tri os for the whole 
population. 



numbers) but also shows that the velocity distribution becomes 
more radially biased with increasing Galactocentric radii. This re- 
sult is in agreement with the findings of cos mological simulations 
(e.g. iDiemand etai]|2007l ; ISales et al.ll2007h . Stars populating the 
outer halo are likely remnants of merger events and therefore can be 
on eccentric orbits. However, the anisotropy parameter /3 is depen- 
dent on the adopted density profile (whereas the rotation parameter 
r\ is independent). More centrally located populations (larger a) ex- 
hibit narrower velocity dispersions (see Section|2j. To compensate 
for the narrower velocity dispersion the unknown tangential distri- 
bution shrinks as we only have a measurement for the line of sight 
velocity. Hence, the distribution becomes more radially anisotropic 
(/3 increases) when we adopt steeper density power laws. This is 
shown in Fig. [5] for two radial bins. To break this degeneracy we 
require a full analysis of the density profile of the BHB sample. As 
our main result regarding rotation is unaffected we defer this task 
to a future paper. 

Evidence for substructure in the metal-rich sample is shown 
in Fig. [6] Here, we show the line of sight velocity dispersion, <7i os , 
as a function of Galactocentric radii. We formulate error bars by 
assuming Gaussian distributed errors and apply a maximum likeli- 
hood routine to find the lcr deviations. There is an obvious 'cold' 
structure in the metal-rich sample for 20 < r/kpc < 30. This 
could be evidence for a shell type structure caused by a relatively 
recent accretion event. The lack of such a feature in the metal-poor 
sample suggests that the net retrograde signal is not a consequence 
of a single accretion event. The dichotomous nature of the outer 
halo may reflect two populations with very different origins. This 
is discussed further in Section|4] 

We note that a larger sample of BHB stars extending to larger 
Galactocentric distances, will be invaluable in studying the outer 
reaches of the stellar halo. To date very few tracers are available, 



but future spectroscopic surveys are likely to increase the sample of 
distant BHBs with measured velocities. For example the proposed 
LSST mission, will provide a vast number of spectroscopic targets 
out to a magnitude limit of r ~ 27. 



3.2 Milky Way globular clusters 

The kinematics of the Milky Wa y globular clusters have been stud - 
ied by a number of authors (e.g. lFrenk & Whitdl 1 98d : IZinrJI 1 985h ■ 
We apply the methods adopted in the previous se ction to mode l the 
halo globular clusters. Using the data taken from lHarrisl d 19961) . we 
restrict attention to those at galactocentric distances r > 10 kpc. 

With line of sight velocities alone, our sample of 41 globular 
clusters poorly constrains /3 and 77. We i mprove this analysis u sing 
the proper motion measure ments from iDinescu et al I dl999h and 
ICasetti-Dinescu et all d2007t) . which include 1 5 of our halo clusters 
(see T able [2]>. We use a similar procedure to IWilkinson & Evansl 
d 19991) and convolve our probabilities with an error function. 
Eqn becomes 



F(l,b,d,vi oa ) = 



£i(vz)£iK) 



x F(l,b,d,vi,v b ,vi os )dvidv b , 



(10) 



where E\ (v) is the error function. We assume the Lorentzian given 
by 



Ei(v) = 



2a( 



s/2kcj\ 20-2 + ^ _ Uoba )2 



(11) 



where u\ is related to the published error estimate by o\ = 
0.477o- meas - The pro perties of such functions a re discussed fur- 
ther in Appendix A of IWilkinson & Ev ans ( 1999). We show in the 
left hand panel of Fig. UJ that the halo clusters have a mildly ra- 



dial distribution (/3 = 0.5 jl '3 ) and net prograde rotation (r\ = 
0.3lo g 5 )- We find a mean streaming motion of (v^) = (58 ± 
27) (r/ 10kpc) _1/4 kms _1 which evaluates to (i^) ~ (60 - 
50) kms -1 in th e range r = (10 — 20) kpc. This a grees with previ- 
ous authors (e.g. lFrenk & White! 1980l : IZinnl 19851) who have found 
that the halo globular clusters are primarily pressure supported, but 
show weak rotation with v TOt ~ 50 — 60 kms -1 . 

The right panel of Fig. |7j shows the v$ velocity components 
of the globular clusters overplotted on a mildly radial model for 
(vcj>(r)). The majority of the globular clusters have positive v,/, and 
hence the net streaming motion is prograde. However, 5 of the halo 
globular clusters have negative v^. As our sample only consists 
of 15 objects, we cannot conclude that this net prograde rotation 
is ubiquitous for the whole population. There is no obvious cor- 
relation between rotational velocity and metallicity, suggesting in- 
dependent accretion events rather than dissipative collapse as the 
formation mechanism. 

The majority of the halo globular clusters in our sample have 
metallicities between —2 < [Fc/H] < — 1 and Galactocentric radii 
in the range 10 < r /kpc < 20 (with the notable exception of Pal 
3). It is interesting that these properties are in common with the 
relatively metal-rich halo BHB stars, which also exhibit net pro- 
grade rotation (albeit a weaker signal). We suggest that some of the 
globular clusters with positive may share a common accretion 
history with these relatively metal-rich field halo stars. 
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Name 


r fknc i 


[Fe/H] 


i'h (kms ) 


fia, cos(5) (mas yr ) 


nx (mas vr ~) 


va, fkms 'l 


Ref 


NGC 288 


11.8 


-1.24 


-46.6 


4.40 ±0.23 


-5.62 ±0.23 


-27 ± 18 


1,2 


NGC 1851 


17.2 


-1.26 


320.9 


1.28 ±0.68 


2.39 ±0.65 


134 ± 29 


1.2 


NGC 1904 


18.9 


-1.54 


207.5 


2.12 ±0.64 


-0.02 ±0.64 


83 ± 29 


1.2 


NGC 2298 


16.0 


-1.85 


148.9 


4.05 ± 1.00 


-1.72 ±0.98 


-27 ± 30 


1.2 


NGC 2808 


11.2 


-1.37 


93.6 


0.58 ±0.45 


2.06 ±0.46 


82 ± 16 


1.3 


NGC 4147 


21.1 


-1.83 


183.2 


-1.85 ±0.82 


-1.30 ±0.82 


67 ±65 


1,2 


NGC 4590 


10.3 


-2.06 


-95.2 


-3.76 ±0.66 


1.79 ± 0.62 


294 ± 30 


1,2 


NGC 5024 


19.0 


-2.07 


-79.1 


0.50 ±1.00 


-0.10 ± 1.00 


240 ± 85 


1,2 


NGC 5272 


12.1 


-1.57 


-148.6 


-1.10 ±0.51 


-2.30 ±0.54 


105 ± 24 


1.2 


NGC 5466 


17.0 


-2.22 


107.7 


-4.65 ±0.82 


0.80 ± 0.82 


-63 ± 64 


1,2 


NGC 6934 


12.4 


-1.54 


-411.4 


1.20 ±1.00 


-5.10 ± 1.00 


-67 ± 60 


1.2 


NGC 7078 


10.5 


-2.22 


-107.5 


-0.95 ±0.51 


-5.63 ±0.50 


129 ± 25 


1.2 


NGC 7089 


10.4 


-1.62 


-5.3 


5.90 ±0.86 


-4.95 ±0.86 


-84 ± 41 


1,2 


Pal 3 


92.9 


-1.66 


83.4 


0.33 ± 0.23 


0.30 ±0.31 


146 ± 95 


1.2 


Pal 5 


17.8 


-1.38 


-55.0 


-1.78 ±0.17 


-2.32 ±0.23 


42 ±34 


1.2 



Table 2. The Milky Way halo globular cluster sample with available proper motions. We give the Galactocentric radii r (note we correct the Harris ( 1996) 
values to a solar position of Rq = 8.5 kpc), the metallicity [F e/H], the helio cent ric velocity, v^, the h elio centric rest frame proper moti ons in right ascension 
and declination and the rotational velocity, v^. References: (l) lHarrisl jl99rj) . (2) lDinescu et alj )l999h . (3)|Casetti-Dinescu et al. 1 2007b. 
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Figure 7. Left: The confidence contours for Milky Way halo globular clusters in the t], /3 plane evaluated from a 2D Maximum likelihood analysis. Arrows 
illustrate prograde, retrograde, radial and tangential distributions respectively. We use — ln(l — /3) on the y-axis to make the ranges occupied by radial and 
tangential models symmetric. The Milky Way globular clusters with available proper motions (N = 15) show evidence for prograde rotation (tj ~ 0.3) and a 
mildly radial velocity distribution (J3 ~ 0.5). Right: The streaming motion as a function of Galactocentric radius for a mildly radial distribution = 0.5 (see 
Fig.[2). Different linestyles represent different rotation parameters (r;). Overplotted are the vj, velocity components of the Milky Way globular clusters with 
available proper motions. The points are colour coded according to metallicity, black being the most metal-poor ([Fc/H] ~ —2.2) and red being the most 
metal-rich ([Fc/H] ~ —1.2). There is no obvious correlation with metallicity but the two clusters with particularly large positive (~ 200 — 300 kms -1 ), 
NGC 4590 and NGC 5024, are also metal poor ([Fc/H] < —2). Note that ~ 5 of the globular clusters have negative v^, but the majority have positive v^. 
As there are only 15 objects in the sample the net prograde motion may not be representative of the whole population. Pal 3 is not shown as its distance of 
r ~ 90 kpc is beyond the range of the plot, but note that it too has a positive (albeit with large formal uncertainties, 146 ± 95 kms ) 



3.3 Milky Way Satellites 

Our final application focusing on our own Milky Way Galaxy looks 
at the satellite galaxies. Fig. [8] shows the result of the maximum 
likelihood analysis keeping the velocity anisotropy fixed at /3 = 0. 
We have split the satellites into classical (L ~ 10 6 Lq) and ultra- 
faint (L ~ 10 2 - 10 4 L Q ) samples (see Tables g] and [3). There 
is some controversy as to whether these groups are separate pop- 
ulations or just the bright and faint components of a single pop- 
ulation. We can see from Fig. [8] that the rotation is poorly con- 



strained. The la confidence limit (horizontal dotted line) encom- 
passes non-rotating, retrograde and prograde solutions. We illus- 
trate the isotropic case (/? = 0) in Fig.[8]but the same form is seen 
for radial or tangential distributions. 

To make further progress, we incorporate the available proper 
motion data of the Milky Way satellites into the analysis. This re- 
stricts the analysis to the classical satellites, as no proper motion 
measurements are available for the ultra-faints. Fig. [9] (left panel) 
shows the resulting likelihood distribution with respect to the rota- 
tion parameter, rj and the anisotropy parameter, /3. 
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Name 


in 


6(°) 


d (kpc) 


«h ( kms ) 


fi a (mas cent ) 


fig (mas cent 1 ) 


vj, (kms 1 ) 


Ref 


Carina 


260.0 


-22.2 


101.0 


224.0 


22 ± 9 


15 ± 9 


68 ± 43 


1,2 


Fomax 


237.1 


-65.7 


138.0 


53.0 


48 ±5 


-36 ±4 


-76 ± 27 


1,3 


Draco 


86.4 


34.7 


82.0 


-293.0 


19 ± 13 


-3 ± 12 


17 ± 47 


1,4 


Leo I 


226.0 


49.1 


250.0 


286.0 








1 


Leon 


220.2 


67.2 


205.0 


76.0 








1 


LMC 


280.5 


-32.9 


50.0 


278.0 


196 ±4 


44 ±4 


82 ±8 


5.6 


Sagittarius 


5.6 


-14.1 


24.0 


140.0 


-283 ± 20 


-133 ± 20 


49 ±22 


1.7 


Sculptor 


287.5 


-83.2 


79.0 


108.0 


9 ± 13 


2± 13 


126 ± 49 


1.8 


Sextans 


243.5 


42.3 


86.0 


227.0 


-26 ± 41 


10 ±44 


193 ± 153 


1.9 


SMC 


302.8 


-44.3 


60.0 


158.0 


75 ±6 


-125 ± 6 


44 ± 16 


5.6 


Ursa Min 


105.0 


44.8 


66.0 


-248.0 


-50 ± 17 


22 ± 16 


-87 ± 48 


1,10 



Table 4. The properties of the classical Milky Way satellites. Here, we give the Galactic coordinates (£, b), the heliocentric distance d, the heliocen tric 
velocity v^, and t he h eliocentric rest fra me proper motions in ri ght ascension and declination an d the rotational veloc ity, v^. References: (l)lMateol 19981 (2) 
Piatek et alj|2003l (3)|piaiek et al Il2007t (4)|Piatek et alj|2008 j . (5) lKarachentsev et al.ll2004 (6) |Piatek et aT]|2008bL (7) lDinescu et alj|2005l (8) |Piatek et alj 
20061 (9) IWalker et alj2008l (10) |Piatek et alj2005l ~ 
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Figure 9. Left: The 2D r\, /3 confidence contour plane for the Milky Way satellite galaxies from the maximum likelihood analysis. Arrows illustrate the regions 
of prograde, retrograde, radial and tangential distributions respectively. The Milky Way satellites with available proper motions show evidence for prograde 
rotation (77 ~ 0.5) and a tangential velocity distribution (/3 ~ —1.5). The inset of the figure shows the confidence contours when the normal vector of the 
angular momentum is tilted with respect to the Milky Way di sc (by an angle 6). More pronounced prograde rotation is seen when 9 » —50°. The dotted line 
shows the approximate tilt angle found bv lMetz et all 120071) using spatial coordinates of the satellites assuming they form a 'disc of satellites'. Right: The 
streaming motion as a function of Galactocentric radius for a tangential distribution fi = —1.5 (see Fig. [2}- Different linestyles represent different rotation 
parameters Overplotted are the velocity components of the Milky Way satellites with available proper motions. Most exhibit positive values but 
Fornax and Ursa Minor have negative . 



With these extra constraints on the tangential velocity compo- 
nents, the maximum likelihood method gives evidence for prograde 
rotation (77 = 0.51q'2> P = — 1.515'q)- Evaluating eqn $6^ gives 
(v^) = (69±42)(r/10kpc) _1/4 kms _1 . Avera ging over the vol- 
ume interval (10 < r/kpc < 150) gives (v^) = (38±23) kms -1 . 
The right hand panel of Fig. [9] shows for the satellites with 
proper motion data overplotted on a tangential model for (v^(r)) 
(see Fig. [2}- Most of the satellites have prograde rotation compo- 
nents, although the Fornax and Ursa Minor satellites show retro- 
grade rotation. 



Wat kins et al.U2010l) recently reported = 0. 44 for the Milky 
Way satellites based on results of simulations by iDiemand et al.1 
l2007t However, based on the proper motion data, the same au- 
thors find P ~ —4.5, favouring tangential or bits. Our result also 
favour s tangential orbits (in agreement with I Wilkinson & Evansl 
dl999h who find/3 w -1). 

The orbits of the satellites for which we have all velocity com- 
ponents are preferentially polar (vg is the la rgest velocity compo- 
nent) . This has been noted in earlier work (e.g lZaritskv & Gonzalez! 
1 19991) and has led to the suggestion that the Milky Way satellites 
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Name 


*( ) 


b ( ) 


d (kpc) i>h (kms x ) 


Ref 


Bool 


358.1 


69.6 


62.0 


99.0 


1 


Boon 


353.7 


68.9 


At* n 


inn 


2 


Can Ven I 


74.3 


79.8 


224.0 


30.9 


3 


Can Ven II 


1 1 3.0 


01 "7 

52. / 


151.0 


-128.0 


2 

J 


Coma 


241.9 


83.6 


44.0 


98.1 


3 


Hercules 


28.7 


36.9 


138.0 


45.0 


3 


Leo IV 


265.4 


56.5 


158.0 


132.3 


3 


LeoV 


261.9 


58.5 


180.0 


173.3 


4 


Segue I 


220.5 


50.4 


23.0 


206.0 


5 


Segue II 


149.4 


-38.1 


35.0 


-39.2 


6 


Ursa Maj II 


152.5 


37.4 


32.0 


-116.5 


3 


Willman I 


158.6 


56.8 


38.0 


-12.3 


1 



Table 3. The properties of the ultra-faint Milky Way satellites. We give the 
Galactic coordinates (I, b), th e heliocentric dista nce d, and the heliocen tric 
velocity v b . Referen ces: (1) iMartin et alJl2007l (2) iKoch et al.| | 2009l (3) 
Simon & Geh3l2007l (4) iBelokurov et"al]|200g. (5) beha et alj|2009l (6) 
Belo kurov et al.l2009l 
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Figure 8. Maximum likelihood analysis of the rotation parameter, r\, for 
Milky Way satellites. The velocity anisotropy is kept fixed at f} = 0. The y- 
axis, Exp[L(r;) — L max ], is the reduction of the likelihood for a particular 
7] value with respect to the maximum likelihood. The red (solid) and blue 
(dashed) lines represent the classical and ultra-faint satellites respectively. 
Rotation is poorly constrained in both cases. 



may occupy a highly inclined disc o f satellites jPalma et~al ] |2002l ; 
iKroupa et ai]|2005l ; (rvletz et al.ll2007h . We have chosen the L z an- 
gular momentum component of the ensemble of satellites to be 
aligned with that of the disc of the Milky Way. As an inset in the left 
hand panel of Fig. [9] we show how the apparent rotation changes 
if we introduce a tilted angular momentum vector for the ensemble 
of satellite galaxies (-90° < 6 < 90°). Note that 9 = cor- 
responds to alignment with the L z angular momentum vector of 
the Milky Way disc. We find more pronounced (prograde) rotation 
when 6 = — 50°l 2 o- This result is not surprising considering we 
find preferentially polar orbits, but does not necessarily mean the 
satellites occupy a rotationally supported d i sc. We also show the 
approximate tilt angle found bv lMetz et ail d2007l) of 8 ss -78° 
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Figure 10. As Fig.|4]but we show the effects of varying the local standard 
of rest velocity (0o). Metal poor ([Fe/H] < -2, N = 1135) and metal 
'rich' ([Fe/H] > —2, N = 2125) subsamples are shown by blue and 
red contours respectively. The paler shades show an increase in the circular 
speed to ©o = 250 kms - 1 whilst the darker shades show a decrease in 
the circular speed to 0o = 190 kms -1 

(dotted line|f| . The authors find the spatial distribution of the satel- 
lites is best described by a highly inclined disc. Whilst our results 
favour a slightly less inclined plane, it is interesting that the de- 
ductions from both a kinematic and spatial analysis broadly agree. 
However, without more accurate proper motion measurements the 
disc of satellites hypothesis cannot be rigorously tested. 



4 THE LOCAL STANDARD OF REST 

Our interpretation of any rotation signal in the Milky Way stel- 
lar halo depends on our assumed local standard of rest rotational 
velocity (Oo). The IAU (International Astronomical Union) rec- 
ommend a value of Oo = 220 kms -1 . In practice, estimates 
of the rotation speed in t he literature vary between 184 kms - 1 
tolling & Merrifieldl 1998b and 272 kms" 1 dMendez et aHll999h . 
Many of these estimates are confined to regions within the solar 
neighbourhood and rely on assuming a value for the Galactocentric 
distance (Ro). 

An abundance of evidence (net rotation, extent, metallic- 
ity, velocity, anisotropy, dispersion profile) suggests the relatively 
metal-rich BHB population of the Milky Way halo is associated 
with the accretion of a massive satellite. The low metallicity and 
flat velocity dispersion profile of the metal-poor population ar- 
gue against an association with a single accretion event. We infer 
that this metal-poor population reflects the 'primordial' metal-poor 
halo. The net retrograde rotation may simply be the result of an un- 
derestimate of the local standard of rest rotational velocity. In Fig. 
[Tol we show the effect of varying Oo on our results for the Milky 



2 There is a typo in the abstract of lMetz et alj fe007l) where an inclination 
angle of 88° is given instead of the (correct) value of 78° 
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Way halo BHB population (see Fig. [4}. Increasing our adopted 0o 
value removes any retrograde signal in the metal-poor population 
and enhances the net prograde signal for the metal-rich population. 
We can provide functional forms for the rotation signal of the two 
components by fitting a simple linear relation to the adopted local 
standard of rest value: 





^275 


"e " 

_220_ 


- 300^ 


kms 


(^t/>)rich = 


^330 


"e ~ 

_220_ 


-315j 


kms - 



(12) 



190 < e /kms -1 < 250 

Here we have removed the dependence on Galactocentric radii 
of (v^) (see eqn[6) by integrating over the appropriate volume in- 
terval (10 < r/kpc < 50). If we assume the metal-poor halo 
component has no net rotation, then we can infer an indepen- 
dent estimate for the local standard of rest rotational velocity of 
6o « 240 kms" 1 . 

An upward revision of the commonly adopted circular s peed 
of 220 kms - agrees with recent work by |Reidetal1 ( l2009h and 
iBovvet all 120091) based on the kinematics of masers found in mas- 
sive star-forming regions of the Milky Way. iReid et al.l J2009I) fit 
for i?o and Qq simultaneously using kinematic and spatial infor- 
mation on sources well beyond the local solar neighbo urhood and 
establi shed a circular rotation speed, Qo ~ 250kms -1 . lBovv et al.l 
re-analysed the maser tracers using less restrictive models 
and found a rotation speed of Oo ~ 240 — 250 kms -1 . They con- 
firm that there is no conflict between recent determinations of the 
circular rotation speed based on di fferent methods (e.g. orbital fit- 
ting of the GD-1 stellar stream, see lKoposov et alj|2010h and quote 
a combined es timate of Qp ~ 240 kms -1 . 

Note that lCarollo et alj J2007I) argued that the net retrograde 
signal i n the outer halo (r > 10 k pc) is due to dynamical friction 
effects. lOuinn & GoodmanHl98 show, based on the effects of dy- 
namical friction, that fragments on retrograde orbits surrender their 
orbital energy to a much smaller extent than those on prograde or- 
bits. The disintegration of many of these fragments by tidal forces 
may result in the formation of a halo stellar popu lation with a net 
retrograde asymmetry (e.g. lNorris&Rvar]|l989h . Hence, in this 
picture the metal-poor component is the accumulation of several 
(smaller) accretion events. However, it is questionable whether this 
is a plausible explanation seeing as dynamical friction mainly af- 
fects massive satellites close to the plane of the disc. This is difficult 
to reconcile with the aggregation of less massive satellites which 
leave debris over a wide range of distances (10 < r/kpc < 50) 
studied in this work. 

We conclude that the lack of a plausible physical explanation 
for a net retrograde rotating metal-poor halo population suggests 
that the commonly adopted local standard of rest circular velocity 
needs to be revised upwards. Assuming the primordial metal-poor 
halo has no net rotation provides an independent estimate for this 
fundamental Galactic parameter. 



5 APPLICATIONS: THE ANDROMEDA GALAXY 

We now apply our analysis to the Andromeda (M31) galaxy. In 
this case, we transform to a coordinate system centred on M31. We 
define the projected distances (in angular units) along the major and 



Name 


Type 


in 


6(°) 


d (kpc) 


i)h kms 1 


Ref 


M31 


Spiral 


121.2 


-21.6 


7851*! 


-301 


1 


And I* 


dSph 


121.7 


-24.8 




-380 


1 


AndH* 


dSph 


128.9 


-29.2 


652li* 


-188 


1 


And IIP 


dSph 


119.4 


-26.3 


7491^ 


-355 


1 


AndV 


dSph 


126.2 


-15.1 


774±g 


-403 


1 


And VI* 


dSph 


106.0 


-36.3 


7831^ 


-354 


1 


And VII 


dSph 


109.5 


-9.9 


763!*! 


-307 


1 


And IX* 


dSph 


123.2 


-19.7 


765l5 5 o 


-207.7 


2 


AndX* 


dSph 


125.8 


-18.0 


702±™ 


-163.8 


3.4 


And XI* 


dSph 


121.7 


-29.0 


760li° 50 


-419.6 


2 


AndXn* 


dSph 


122.0 


-29.5 


8301 1 ™ 


-558.4 


2 


And XIII* 


dSph 


123.0 


-29.9 


910+ 30 


-195.0 


2 


AndXIV* 


dSph 


123.0 


-33.2 


740 +no 


-481.1 


5 


And XV* 


dSph 


127.9 


-24.5 


770+™ 


-339 


6 


And XVI 


dSph 


124.9 


-30.5 


5251!° 


-385 


6 


Pisces 


dlrr/dSph 


126.8 


-40.9 


76911 


-286 


1 


Pegasus 


dlrr/dSph 


94.8 


-43.6 


oiq+30 

a -30 


-182 


1 


NGC 147 


dE 


119.8 


-14.3 


6751?, 7 . 


-193 




IC 1613 


dE 


129.8 


-60.6 


7001^ 


-232 




NGC 185 


dE 


120.79 


-14.5 


616l*« 


-202 




NGC 205 


dE 


120.7 


-21.1 


8241* 7 , 


-244 




IC 10 


dhT 


119.0 


-3.3 


8251!° 


-344 




M32 


cE 


121.2 


-22.0 


7851*! 


-205 





Table S. Properties of the Andromeda Satellites: Here, we give the satel- 
lite name and type (dSph = dwarf spheroidal, dlrr = dwarf irregular, dE = 
dwarf elliptical and cE = classical elliptical), Galactic coordinates (I, b), he- 
liocentric distance d and heliocentric velocit y v^. Satellites belonging to 
the rotating group are s tarr ed. References: (l)lMcCon nachie & Irwinl2o"o6l 
(2)ICollins etal.ll2009l OllZucker & et ajl2007l (4) lKalirai et alj|2009l. (5) 
iMaiewski et al.l2007l (6) lLetarte et alj2009l 



minor axis, X and Y via: 

xo = sin(a — Qo)cos(c5), 
yo = sin(<5)cos(<5o) — cos(a — ao)cos(<5)sin((5o), 
X = x sin(# a ) + y cos(6 l a ), 
Y = -x cos(6» a ) + y sin(0 a ), 

where xo,yo are Cartesian coordinates for a, S with respect to 
the centre of M3 1 (a = 00 h 42 m 44 3 s So = +41°16 / 9"; 
ICotton etal.lll999h and 6 a = 37.7° dde Vaucouleursl 1 19581) is 
the position angle. We adopt the tran sformation outlined in Ap- 
pendix A of Evans &deZeeuwl ( ll994l) to relate the projected co- 
ordinates to the galaxy coordinates (assuming an inclination angle, 
% = 77.5°). 

5.1 M31 Satellites 

The main drawback in the M3 1 analysis is that the distances to the 
satellites are not as well constrained as for the Milky Way satel- 
lites. To take this into account, we integrate the distribution func- 
tion along the line of sight between the distance errors. The left 
hand panel of Fig. QT] shows the results of the maximum likeli- 
hood analysis applied to all of the M31 satellites (see Table |3)- 
Unlike the situation for the Milky Way satellites, our line of sight 
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to M3 1 is not constrained to a predominantly radial velocity com- 
ponent. Fig. [TT] shows evidence of prograde rotation from the line 
of sight velocity alone. We further investigate this rotating com- 
ponent by splitting the M31 sample into two groups based on the 
satellite type. These are shown by the red (dSph) and blue points 
(non-dSph) in Fig. II II where we find a rotating group (dSph, red) 
and a non-rotating group (non-dSph, blue). The rotating subset 
has a rotation parameter, r\ — 0.4lgJ and anisotropy parameter, 
[3 — O.llg'5, which corresponds to a streaming motion, (v^) — 
(62 ± 34)(r/10kpc)" 1/4 kms" 1 . Averaging over a suitable vol- 
ume interval (10 < r/kpc < 200) gives (v^) = (32± 17) kms~\ 

The observed line of sight velocities of the M31 satellites can 
be decomposed into radial and tangential components by making 
use of the spatial information relative to M31. These are lower 
limits as we only have the line of sight velocity component (see 
lMetzetal.ll2007h . There are also large uncertainties on the dis- 
tances of the satellites which permeate into the spatial coordinates 
with respect to M3 1 . However, we still find it instructive to look at 
these approximate velocity components. The right panel of Fig.fTTI 
shows the lower limit tangential and radial velocity components 
for the M31 satellites. The members of the dSph rotating group are 
shown by red stars and the non-rotating satellites are shown by the 
blue circles. The majority of the rotating subgroup have a larger 
tangential velocity than radial velocity component. 

The rotating subgroup consists entirely of dwarf spheroidal 
galaxies. Satellites with similar accretion histories can share com- 
mon properties such as their 'type' (i.e dSph) and orbital orienta- 
tion. The apparent rotation of the dSph satellites of M3 1 suggests 
that these dwarfs may share a common origin. 



5.2 M31 globular clusters 

Over recent years, studi e s by several aut h ors (e.g. IPerrett et al.l 
120021 iGalleti et alj 12001 iKim et all 12001 ICaldwell et al.l |2009|) 
have produced a rich sample of M3 1 globular clusters with avail- 
able kinem atic data. Our data derives from the Revised Bologna 
Catalogue jGalleti et all |2004) . Distance measurements are un- 
available for the globular clusters, so we integrate our distribution 
function along the line of sight. 

We consider the 93 globular clusters for which the projected 
distance along the semi-minor axis is gr eater than 5 kpc (\Y\ > 
5 kpc). This criterion was introduced bv lLeeetalJ ( 120081) . For 54 
of these clusters, metallicities are available and so we further di- 
vide the sample into 24 metal-rich ([Fe/H] > —1) and 30 metal- 
poor ([Fe/H] < — 1) components. Fig. 1 121 shows the outer M31 
globular clusters are rotating prograde with respect to the M31 
disc. The metal-rich subset (solid red contours) has r\ — OIq'j 
and ft = 0.3155'g gi ym g a streaming motion, (u^) = (94 ± 
35)(r/10kpc)- 1/4 kms- 1 . The metal -poor subset (dashed blue 
contours) has a less pronounced rotation with r\ = 0.1 ± 0.5, 



and a more pronounced radial anisotropy, ft = 0.71q'2- This in- 
dicates that it has a smaller streaming velocity: (v^,) — (53 ± 
34)(r/10 kpc)~ 1//4 kms -1 . However, the broad overlap of the lcr 
error envelopes of these samples mean that it is premature to con- 
clude that we are indeed seeing a signifi cant diffe re nce be tween the 
metal-rich and metal-poor populations. iLee et alj d2008l) also find 
strong rotation for the outer globular clusters, suggesting a hot ro- 
tating halo or an extended bulge system in M3 1 . The role of the halo 
in M31 has been contentious. Some stu dies have suggested M31 




Figure 12. The confidence contours of the orbital parameters 77 and /3 for 
the M31 globular clusters. We only consider globular clusters where \ Y | > 
5 kpc to probe the outer regions and reduce disc contamination. The solid 
red and dashed blue contours represent the metal-rich ([Fe/H] > -1) and 
metal-poor ([Fe/H] < -1) subsamples respectively 



has a n extended bulge but no ha l o (e.g. IPritchet & van den Bergl 
1 1994 iHurlev-Keller et alj |2004 llrwin et alj 120051 : iMerrett et al 



I2006I) , whilst more recent studies have revealed evidence for a huge 
halo JChapman et"aT]|2006J : iKalirai et alJlioOol : llbata et al.1 120071) . 
The weaker rotational component in the metal-poor sample may 
suggest a transition from the bulge to a halo component but a 
larger kinematic sample extending further along the minor-axis is 
required for confirmation. 



6 CONCLUSIONS 

We studied the rotational properties of halo populations in the 
Milky Way Galaxy and Andromeda (M31). We modelled the grav- 
itational potential, $, and the density profile of the tracer popula- 
tion, p, as spherically symmetric power-laws with indices 7 and a 
respectively. The shape of the velocity distribution is controlled by 
two parameters, namely f3 which governs the velocity anisotropy 
and r\ which governs the rotation. 

For our applications, we used the potential $ ~ r 1 with 

7 = 0.5 as this is a good approximation to a Navarro-Frenk- White 
profile at large radii. Our tracer populations have a density profile 
p ~ r~ a with typically a — 3.5. A maximum likelihood method 
is used to constrain the orbital parameters /3 and rj for five differ- 
ent populations: (1) Milky way halo blue horizontal branch (BHB) 
stars, (2) Milky Way globular clusters, (3) Milky Way satellites, (4) 
M31 satellites and (5) M31 globular clusters. We summarise our 
conclusions as follows: 

(1) From our constructed sample of halo BHB stars from SDSS, 
we find evidence for two different halo populations: a pro- 
grade, rotating, comparatively metal-rich component ([Fe/H] > 
—2) and a retrograde, rotating, comparatively metal-poor com- 
ponent ([Fe/H] < —2). The mean streaming motions are 
(v*) = (21 ± ll)(r/10kpc)~ 1/4 kms- 1 and (vj,) = -(35 ± 
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V 

Figure 11. Left: The confidence contours of the orbital parameters r\ and /3 for the M31 satellites. The red and blue points correspond to the dSphs (rotating) 
and non-dSph (non-rotating) subsets respectively. Arrows represent the regions of prograde, retrograde, radial and tangential distributions. Right: The line of 
sight velocities of the M3 1 satellites decomposed into tangential and radial components. These are lower limits as only one component of velocity is known. 
The red stars and blue circles represent the members of the dSph and non-dSph groups respectively. 



10) (r/10 kpc)~ 1,/4 kms -1 for the richer and poorer populations 
respectively. A rough estimate suggests that these structures may 
have stellar masses of M ~ 1O 8 M0. We find evidence from the 
velocity dispersion profile that the metal-rich sample is associated 
with the accretion of a massive satellite. The line of sight velocity 
dispersion of the metal-rich sample has a kinematically cold com- 
ponent at 20 < r/kpc < 30, which could be evidence for a shell 
like structure caused by the accretion event. We find no such evi- 
dence to suggest the metal-poor sample is associated with an accre- 
tion event. We suggest that this population typifies the primordial 
stellar halo and the net retrograde rotation reflects an underestimate 
in our adopted local standard of rest circular velocity, which may 
be as high as ~ 240 kms -1 . 

(2) The halo globular clusters of the Milky Way show evidence of 
net prograde rotation, (u^) = (58 ± 27)(r/10 kpc) -1 ' 4 kms -1 , 
and have a mildly radial velocity distribution. Previous studies have 
found a mild net prograde rotation from an analys is of the line of 
sight velocities of the halo globular clusters (e.g. iFrenk & White! 
Il980h . Our results are in agreement, but this deduction arises from 
the analysis of those 15 of the 41 halo globular clusters which have 
available proper motions. We note that the majority of our globular 
cluster sample have metallicities in the range —2 < [Fe/H] < — 1 
and are located between 10 < r/kpc < 20. It is interesting that we 
also see net prograde rotation in the halo BHB stars in this distance 
and metallicity range. This suggests the relatively metal-rich BHBs 
and halo globular clusters may share a similar formation history. 

(3) Line of sight velocities alone poorly constrain the orbital param- 
eters of the Milky Way satellites. Using available proper motions 
for 9 of the classical satellites, we find evidence for prograde ro- 
tation, (v^) = (69 ± 42)(r/10kpc) _1/4 kms -1 , and a tangential 
anisotropy parameter /? = — 1.5j^o. The orbits are preferentially 



polar. This leads to the detection of a more pronounced rotation 
when the normal vector of rotation is inclined by \8\ « 50° to the 
normal of the disc. This hints that the satellites may be part of a 
rotationally supported disc (as postulated by a number of authors), 
but without more accurate proper motions this hypothesis cannot 
be confirmed. 

(4) Application of these methods to the M31 satellites is po- 
tentially powerful, as the line of sight velocity has contributions 
from both the tangential and radial velocity components with re- 
spect to the M31 centre. The overall population shows hints of 
prograde rotation. We can divide the satellites into two groups. 
One possesses pronounced prograde rotation, (v^) = (62 ± 
34)(r/10 kpc)~ 1,/4 kms -1 , whilst the other shows no evidence for 
a rotating signal. Interestingly, the rotating group consists entirely 
of the dwarf spheroidal satellites. The correlation between satellite 
type and orbital orientation suggests these satellites may share a 
common origin. 

(5) There is a rich dataset on the M31 globular clusters which 
remains relatively unexploited. We examined the globular clus- 
ters with projected distance along the minor-axis greater than 5 
kpc. These outer globular clusters have prograde rotation which is 
more pronounced for the metal-rich ([Fe/H] > -1) subset. We find 

= (94±35)(r-/10kpc) -1/4 kms -1 for the metal-rich popu- 
lation and (v$) = (53 ± 34)(r/10 kpc) -1/4 kms -1 for the metal- 
poor population. The strong rotation, especially in the metal-rich 
sample, suggests these globular clusters may belong to the M31 
bulge system. However, a more robust explanation awaits a more 
substantial kinematic sample extending further along the minor- 
axis. 

The Gaia satellite promises to provide proper motion data on all 
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Figure 13. The estimated error in the mean streaming motion of the Milky 

Way satellites as a function of number of satellites. We assume random 

eiTors on the tangential velocity components in the range 

15kms -1 . Three different models are used, shown by the blue (rj = 0.9, 

= 0.1), red (77 = 1.9, & = 0.6) and green (77 = 0.5, & = -0.7) lines 

respectively. 



dwarf gal axies of the Milky Way an d M31 and for thousands of 
halo stars. Iwilkinson & Evansl i 19991) suggested that Gaia can con- 
strain the proper motion of distant classical dwarfs in the Milky 
way, like Leo I, to within ±15 kms -1 and nearer dwarfs like Ursa 
Minor to within ±1 kms - . Whilst Gaia will improve the proper 
motion measurements for the classical satellites, deeper comple- 
ments to Gaia such as the LSST will be needed to achieve this 
level of accuracy for the ultra-faint satellites. 

In anticipation of this future work, we close by consider- 
ing Monte Carlo simulations of our distribution functions to es- 
timate how well we can constrain rotation with these future as- 
trometric surveys. We assume random errors on the tangential 
velocity components between lkms -1 and 15 kms -1 . Fig. 1131 
gives a rough estimate of how the accuracy depends on the num- 
ber of satellites with proper motion measurements. For N w 
20 satellites, we can constrain the mean streaming motion to 
cr v w 30(r/10 kpc) _1//4 kms -1 . This evaluates to a v . ~ (10 — 
20) kms -1 for satellites between r = (50 — 250) kpc. The com- 
bination of more accurate proper motion measurements with radial 
velocities can potentially probe any apparent kinematic dichotomy 
between the classical and ultra-faint satellites, and perhaps provide 
more robust evidence bearing on the disc of satellites. 

Our analysis of the Milky Way stellar halo using BHB stars 
will be greatly improved by increasing the number of targets at 
larger radii. Coupled with more accurate radial velocity measure- 
ments, we will be able to map the kinematic structure of the stel- 
lar halo in greater detail and see if the apparent dual nature of 
the metal-rich and metal-poor components continues out to further 
reaches of the halo. LSST will hopefully play a vital role in finding 
these BHB targets. 
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